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SUMMARY 


A generalized methodology to predict the fatigue life and reliability of 
a rotating disk such as used for aircraft engine turbines and compressors Is 
advanced. The approach Incorporates the computed life of elemental stress 
volumes to predict system life and reliability. Disk speed and thermal gra- 
dients as well as design variables such as disk diameter and thickness and bolt 
hole size, number and location are considered. 


INTRODUCTION 

A key to cost effective design Is the ability to accurately predict com- 
ponent life without the need for extensive testing or field experience. 

Current military support programs for gas turbine engines depend on life pre- 
diction methods to plan overhauls and part replacement. Errors In life pre- 
diction result In premature part replacement or Insufficient availability of 
spare parts (ref. 1). 

Life prediction methods for most machine elements are typically based on 
yield stresses and fatigue limiting stresses. However, In addition to the 
component's material properties, proper consideration must be given to the 
effects of notches, surface condition, component size, residual stress, tem- 
perature, duty cycle and environmental factors. Problems can also arise from 
lack of reliable fatigue data, overly-conservatlve design factors. Inconsist- 
encies among computer codes, and misinterpretation of the actual design loads. 
To compound these difficulties, the established fatigue limit for much of the 
reported data Is a mean or average value. This can be Interpreted as meaning 
that life Is finite, even at the fatigue limiting stress (ref. 2). 

In aerospace applications, an engineer must be especially cognizant of 
size and weight constraints which Impact design decisions. While designing at 
or below the material fatigue limit may be desirable In most Industrial appli- 
cations, in aerospace applications. It Is almost mandatory to design for a 
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finite life at an acceptable probability of survival. Zaretsky (ref. 2 ) out- 
lined such a methodology based In part on work of W. Welbull (refs. 3 and 4) 
and G. Lundberg and A. Palmgren (refs. 5 to 7). Zaretsky's method (ref. 2) Is 
similar In approach to that of Ioannldes and Harris (ref. 8 ). 

Mahorter, et al . (ref. 1 ), discuss the accuracy of life prediction tech- 
niques for military turbine engine components such as compressor and turbine 
disks. The development of a 0.794 mm (0.03 In.) crack In any of the critical 
areas of a disk, such as bolt holes, bore, dovetail, etc.. Is considered the 
end of Its low cycle fatigue (LCF) life (ref. 1). This Is the life, Lqj, 
at which the military requires a 99.9 percent probability of survival or a 0.1 
percent failure rate. Disk retirement policy requires that the disks be 
removed from service or reworked at their Lq.i life (ref. 1 ). The statis- 
tically predicted lives for bolt hole cracks did not agree with the analyti- 
cally predicted values. These results Imply that a deterministic approach to 
life predictions Is not necessarily conservative, and that a probabilistic 
approach Is viable In light of all the statistical variations In the design 
parameters. 

In view of the aforementioned. It Is the objective of this work to 
(a) apply the method of Zaretsky (ref. 2) to statistically predict the life of 
a generic solid disk with and without bolt holes, (b) determine the effect of 
disk design variables, thermal loads, and speed on relative life, and 
(c) develop a generalized equation for determining disk life. 


NOMENCLATURE 

A Material constant, stress cycles 

c stress-life exponent 

D disk diameter, m (In.) 

E Material Young's Modulus, N/m^ ( ps 1 ) 

e Welbull slope 

Kq proportionality constant, m (In.) 

Kl proportionality constant, m^/N (psl -1 ) 

Kn proportionality constant, rpm 

K-j- proportionality constant, K (°F) 

Kt proportionality constant, m (In.) 

L life, stress cycles 

L 0 characteristic life, or life at which 63.2 percent of population fall 
L s system life, stress cycles 
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L-| o ten percent life, life at which 90 percent of a population will survive, 
stress cycles 

i 

N speed, rpm 

N 0 Dynamic Speed Capacity, rpm 

r radial coordinate from disk centerline, m (In.) 

R disk radius, m (In.) 

S probability of survival 

t disk thickness, m (In.) 

T steady-state temperature, K (°F) 

V stress volume, m^ (In. 3) 

a coefficient of thermal expansion, m/K (1n./°F) 

<»L fatigue limit stress, N/m 2 (psl) 

o rr radial stress component, N/m 2 (psl) 

Dee tangential stress component, N/m 2 (psl) 

t shear stress, N/m 2 (psl) 

tl fatigue limit stress, N/m 2 (psl) 

Subscripts: 

a,b states a and b, or bodies a and b 
1 ith component, and denotes Inner disk radius 
o reference point, and denotes outer disk radius 

ANALYSIS 
Failure Theory 

The Welbull Weakest Link Theory (refs. 3 and 4) states that for a body, 
the probability of survival, S, Is a function of at least five parameters. 
These five parameters are; fatigue life, L; material Welbull slope, e; stress 
state, a; stress-life exponent, c; and stressed volume, V. In mathematical 
form this relation Is expressed as: 

lnl~o C L e V (1) 
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Successful Implementation of equation (1) requires that three effects be 
addressed. These are: 

(1) stress/llfe characteristics 

(2) statistical failure distributions 

(3) stressed volume variations 

A means In which the stress-life characteristic of a material can be 
defined Is found from the stress-life relation general expression where, 

L = A <T C (2) 


and where A Is a material constant. The values of A and c In this rela- 
tion can be determined from experimental fatigue data obtained from coupon 
tests (ref. 9). 

Rearrangement of equation (2) will allow life prediction at a certain 
stress state given knowledge of a reference state where the stress and cycles 
to failure are known. The equation 



gives life Lb at a known loading as related to state a. This equation 
will allow life estimation given knowledge of the stress-life exponent c and 
a reference state a. (Note that the constant A from eq. (2) Is not required 
In eq. (3).) 

The specification of material fatigue characteristics cannot be solely 
defined through the use of equation (3). The material data obtained for 
determining material values In equations (2) and (3), are usually presented as 
mean or median values. To account for deviations from this mean value, a 
Welbull distribution Is used to describe the failure distribution. The Welbull 
distribution Indicates the statistical probability of survival of a given 
material or structure, under a particular loading, as a function of the number 
of cycles subjected to that loading. The determination of the Welbull distri- 
bution for loading conditions other than those for which test data Is available 
will require knowledge of the Welbull slope and a point on the new Welbull 
line. The equation of a line using Welbull coordinates Is: 



for a single body at constant loading with S 0 and L 0 as known reference 
points. 

An Important Implication of equation (1) Is that as the volume over which 
a critical stress Is acting Is Increased, a reduction In survivability will 
occur for the same number of cycles. Manipulation of equation 1 for two vol- 
umes subjected to equal stress gives: 
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( 5 ) 


S. = S a 

d a 


v h /v a 

b a 


for constant values L, e, and c. This equation shows that for two bodies of 
same material and subjected to Identical loading, the body of larger volume 
will contain more probable defects and, thus, have a lower probability of sur- 
viving a fixed number of stress cycles. Equation (2) has been supported by 
much experimental evidence (ref. 2). 

A way for predicting the life of a general body "a* related to the life 
of body M b" Is given by: 



The derivation for equation (6) Is given Is appendix A. This equation Is sim- 
ilar In form to equation (5), but also Includes the effect of differences In 
stressed volume. 

If n components combine to form a mechanical system, the system relia- 
bility relation states that: 


n 



where L s Is the total system life and Is the 1*h component's life 
(ref. 2). For this formulation, each component may have a different life but 
the survivability Is constant for all the components, and hence for the system. 

Methodology 

The concept outlined In equation (7) can easily be extended from a 
mechanical system of discrete components to a continuous structure. This pro- 
cedure can be used with the results of a finite element analysis, whereby a 
complex structure Is modeled by discrete elements. Each element of the model 
will have a distinct volume and load-dependent stress. These elemental char- 
acteristics can be applied through equation (7) to determine the elements' 
lives and reliabilities and thus, the structure's life and reliability. The 
Implementation of equation (7) to determine structural life and reliability Is 
shown In the analysis methodology flowchart of figure 1. 

The first step to fatigue life determination Is to model the structure 
Into elemental volumes and stresses. The general guidelines for finite element 
modeling can be used for this structure modeling. In this manner, the ele- 
mental stresses and volumes can be used directly within equations (6) and (7). 

The second step In the fatigue life analysis Involves selecting a critical 
element within the finite element model based on the maximum stress state. 
Because of the formulation of equations (1) to (5), all the element stresses 
and volumes are normalized with respect to the critical element. 
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Several different criteria can be considered for normalizing the critical 
structural element. These criteria are: 

(1) Maximum Shear Stress Failure Theory 

(2) Maximum Normal Stress Failure Theory 

(3) Maximum Element Strain Energy Density 

(4) Maximum Element Volume 

The maximum shear stress should generally be used In Identifying the 
critical element. Justification Is based on the applicability of the maximum 
shear stress failure theory for ductile materials. Independent of the selected 
criteria, for all calculation purposes the element with maximum critical stress 
Is given a relative life of unity at a probability of survival of 90 percent, 
or a 10-percent probability of failure. 

The normalized stresses and volumes are used to compute elemental lives 
relative to the critical element. The relative lives are combined using equa- 
tion (7) to give a structural life relative to the critical element. Two key 
advantages of using this relative life methodology are: 

(1) Simplified qualitative comparison between differing designs 

(2) Reduced fatigue life data requirements 

By normalizing a family of candidate designs through a singular critical ele- 
ment, qualitative assessments of designs changes on structural life can be 
easily and quickly made. By specifying relative lives, only one series of 
fatigue coupon tests need be run to establish the necessary Welbull parameters. 


RESULTS AND DISCUSSION 
Solid disk 

Parametric analytical studies were conducted to Investigate the effect of 
varying physical dimensions and speed on the relative lives of a generic solid 
disk. The physical model of the solid disk requires that the structure be 
divided Into a series of components (elements). This Is most easily done by 
considering the disk as a collection of concentric rings. The disk was modeled 
with 10 rings of equal radial Increment for the parametric values of disk 
diameter studied. Expressions for the average hoop, radial, and maximum shear 
stresses occurring across the rings are derived for equations given In 
(ref. 10). 

Assuming a stress-life exponent of 9 (refs. 8 and 9), the effect of 
Welbull slope modulus, on life was determined. The results are shown In fig- 
ure 3. Welbull slopes of 1, 2, and 3.57 are equivalent to exponential, 
Rayleigh, and normal (Gaussian) distributions respectively (ref. 2). For 
Welbull slopes of approximately 1.5 or greater, life Is fairly constant. For 
Welbull slopes between 1 and 1.5 life Increases by approximately 15 percent. 

As a result, for the remainder of the analysis this effect will be considered 
negligible. 

The effect of disk diameter on life for a speed of 9000 rpm and a stress- 
life exponent, c, of 9 Is shown In figure 4. As the disk diameter Is 
Increased, both the stress and the stress volume Increase and life will 
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decrease. The effect of Welbull slope Is negligible. The equation for this 
relation for the effect of diameter on life Is 



where D Is the disk diameter, and Kq Is a constant equal to 0.61 m 
(24.0 In.). 


( 8 ) 


The effect of disk thickness on the Lio life Is shown In figure 5. 
Because the stressed volume Is Increased, life will decrease. The effect of 
Welbull slope Is negligible. The equation for the relation for the effect of 
disk thickness on life Is as follows: 




(9) 


where T Is the disk thickness, and Kt Is a constant equal to 0.0254 m 
(1.0 In.). 

As disk speed Is Increased, stresses within the disk will Increase. These 
stress Increases will cause a decrease In life. The effect of disk rotational 
speed on L^o life is shown In figure 6. The equation for the effect of 
speed on life be written as follows: 


L 


10 



14.3 


( 10 ) 


where N Is the disk speed and Is a speed constant equal to 9000 rpm. 

For each stressed elemental volume within the body of the disk, the life 
was determined using an Inverse stress life relation shown In equations (2) 
and ( 3) C Not all materials will exhibit the same stress-life relation. The 
stress^llfe exponent Is generally determined experimentally. For the previous 
calculations, a stress-life exponent of 9 was assumed. Using varying values 
of c with a reference disk, the effect of the stress-life exponent on the 
L-jq life of the disk was determined. These results are shown In figure 7. 

The equation for the effect of the stress-life exponent on disk life Is as 
follows: 


L 


10 


0 


0.606 


(ID 


A concept first expanded by Ioannldes and Harris (ref. 8) was the appli- 
cation of a fatigue limit In the determination of the lives of the elemental 
stressed volumes. Basically, the concept Is that where the shearing stress Is 
equal to or less than the value determined or assumed for the fatigue limit, 
the probability of survival for that elemental stressed volume Is 100 percent. 
In their examples, Ioannldes and Harris (ref. 8) use the following relation: 


7 


1 


( 12 ) 


where o Is the operating stress and ol I s the fatigue limiting stress at 
or below which no failure Is expected to occur. 

The relation represented by equation (12) will result In a stress-life 
exponent Inconsistent with equations (2) and (3) and was not used the analysis 
reported herein. Instead, equations (2) and (3) were used In determining the 
probability of survival and the resultant life for values of t < A 
100 percent probability of survival was assumed. The results of assuming 
varying values of a fatigue limit for the shearing stress Is shown In figure 8. 
An approximate relation between life and tl can be written as 

L] o ~ exp(KL • tl) (13) 

where Kl Is equal to 4.02x10-8 m 2 /N (2.77x1 0~ 4 ps 1 “1 ) . 

The effect of a radial temperature difference upon the overall disk rel- 
ative life was determined by considering the thermal stresses superimposed upon 
the centrifugal disk stresses due to disk rotation. A steady-state radial 
temperature distribution was applied to a disk with a very small central hole. 

A disk with a small central hole was considered because this configuration 
relieved a mathematical singularity which occurs for Isothermal boundary con- 
ditions upon a solid disk. The solution of the heat diffusion equation for 
steady-state, one-dimensional heat transfer yields a temperature distribution 
as: 


T = 



In 



(14) 


for T^ and T 0 the Inner and outer radii temperatures. Assuming that 
small temperature differences exist (AT < 500 °F), an elasticity solution for 
the disk thermal stresses was used as: 



(15) 


(16) 


where T Is the steady-state change In temperature from the Initial stress- 
free temperature found from equation (14). The total disk stresses due to both 
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the Centrifugal disk loading and the thermal stresses can then be found by 
superposition. 

This method Is outlined In appendix B. The effect of various radial tem- 
perature differences on disk relative life Is shown In figure 9. An expression 
for this thermal effect Is as follows: 


L 


10 





(17) 


where AT Is the radial temperature difference and Kj Is constant equal 
0.56 K where AT Is In Kelvin or IF where AT Is In Fahrenheit and where 
AT Is equal to or greater than Kj. 

By combining the product of equations (8) to (11), (13) and (17), the 
following equation for the 10-percent life of a uniform solid rotating disk may 
be obtained: 




(exp (K.,t l )) 


(18) 


Applying the concept of a material factor. A, which can be determined from 
experiment or from field experience, equation (14) can be rewritten as follows: 


10 = A 




(exp (K-j t l ) ) 


(19) 


Introducing the concept of a Dynamic Speed Capacity N 0 which Is defined 
by Zaretsky (ref. 2) as the speed which would produce a theoretical life of one 
million stress cycles, then from equation (19), 









V, 

AxlO" 6 

W f 

t\ „ 14.3 /9 

7 N V 

^0.606 ^ 

K \0.52 

A?) < ex P ( k l t i }) 



1/14.3 


( 20 ) 


and for any speed N for a given disk geometry, 

14.3 

/V 

U 


-10 


(?) 


( 21 ) 


From the Welbull equations (3 and 4), 


in in (1 ) 


e In 



( 22 ) 


and using values of S equal 0.9 and L equal L-j q from equations (19) or 
(21), the characteristic life Lg, or the life at which 63.2 percent of the 
disks can be expected to fall, can be determined. Having the cha acterlstlc 
life, any life at any given probability of survival, such as the life at a 99.9 
probability of survival, can be determined from equation (22). 
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Effect of Bolt Holes 


As discussed by Mahorter, et al , (ref. 1), failure of the disk will occur 
at tie-bolt holes In the disk. This will be the critical location for failure 
on a disk. Using the finite element model shown In figure 2, the effect on 
disk life of bolt hole size, location, and number were determined. The results 
of this analysis are shown In figures 9 to 11. The relative lives shown were 
normalized relative to the life of a solid disk with no holes having a 0.61 m 
(24.0 In.) diameter and being 0.0254 m (1.0 In.) thick, running at 9000 rpm 
with a Welbull slope or modulus of 1.00 and a stress-life exponent, c, of 9. 

Figure 10 shows the effect of bolt diameter on disk fatigue life. It 
would appear that at bolt holes having a diameter of less than 10.2 mm 
(0.4 In.), the effect of hole size Is nominal. At bolt holes larger In diam- 
eter, the effect Is most significant. This analysis would suggest that the 
bolt holes In a disk should be smaller but more numerous In number. 

The effect of Increasing the number of bolt holes Is shown In figure 11. 
Increasing the number of holes and keeping the hole diameter less than 10.2 mm 
(0.4 In.) appear to have less of an effect on life than having fewer holes 
larger In diameter. Figures 10 and 11 can be used to obtain life modifying 
factors for equations (19) or (21). 

Figure 12 shows the effect of hole location In the disk as a percent of 
the disk radius measured from the axis of rotation. The results show a general 
trend of Increasing life as the holes are moved radially outward. This would 
not be unexpected since stresses decrease with Increasing distance from the 
center. However, the trend Is not totally convincing since at a location 
between 40 and 50 percent of the disk radius, there Is an Inversion of the 
trend, that Is, life decreases and then begins to Increase again. It cannot 
be determined with reasonable certainty that this Is correct or that the finite 
element mesh size was properly selected for these calculations. Further anal- 
ysis Is required. However, It can be reasonably concluded that the bolt holes 
should be placed at a great distance from the center of the disk as Is 
practical . 

As discussed In the previous section, the effect of Welbull slope was 
considered negligible for the solid disk with no holes. However, for the disk 
with bolt holes, there appears to be a significant reduction In life with 
Increases In slope. This trend was not expected. The size of the finite ele- 
ment mesh may have contributed to this phenomenon. However, the effect of mesh 
size was beyond the scope of the present Investigation. 


SUMMARY OF RESULTS 

A method was developed to predict the fatigue life and reliability of a 
rotating disk such as used for aircraft engine turbines and compressors. The 
method Incorporates the computed life at a given probability of survival of 
elemental stressed volumes within the body of the disk to arrive at total disk 
life and reliability. The effect on life of disk speed as well as design var- 
iables such as disk diameter and thickness and bolt hole size, number and 
location were determined. An Investigation Into the effect of a radial tem- 
perature difference and the accompanying Induced thermal stresses on disk rel- 
ative life was also performed. In addition, material characteristics such as 
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Welbull modules and fatigue limit were considered. A Dynamic Speed Capacity 
N 0 was Introduced. The following results were obtained: 

1. The life of a solid disk was Inversely proportional to speed to the 
14.3 power, diameter to the 20 power, thickness to the first power, and radial 
temperature difference to the 0.52 power. 

2. For assumed material characteristics, life was Inversely proportional 
to the stress-life exponent to the 0.606 power, and exponentially related to 
the fatigue limiting stress. 

3. Increasing the number of bolt holes and keeping the hole diameter less 
than 10.2 mm (0.4 In.) appear to have less of an effect on life than having 
fewer holes larger In diameter. 

4. As hole location Is moved out radially, disk life Increases suggesting 
placement of bolt holes at the furthest practical radial location. 
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APPENDIX A 


The equation for a straight line on Welbull paper Is: 


In 



In In — = e( InL. 
^a 


" ln L a > 


which simplifies to: 



(A.l) 


This equation represents the survival probability for body "a" for any desired 
life at a fixed applied stress o a . 

In order to develop an expression for a new load "b", estimates of the 
associated and Lb must be made. These are: 

s b = s a 


and, 


L h = L 
b a 


ft) 


(A. 2) 


Now, If body "b" has stress ob, but a different volume Vb- 


Sr = S, 

b a 


v a /v h 

a b 


(A. 3) 


Using equations (A.l), (A. 2), and (A. 3), and equal element survivabil- 
ities, S = S a , the reference life Lb can be given by: 


i ■ (?)' ft) 


l/e 


(A. 4) 
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APPENDIX B. - Thermal Analysis 


The solution of the steady heat diffusion equation In cylindrical coor- 
dinates will give a steady-state temperature distribution of [11]: 


T(r) 



(B.l ) 


After Isothermal boundary conditions are applied at a critical Inner radius 
, and the disk outer radius R 0 . A critical Inner radius Is necessary 
to relieve a mathematical singularity at the center of a solid disk. The tem- 
perature within the Inner radius was assumed constant at T^. 

An elasticity solution for the thermal stresses Induced within a thin disk 
will require solving the following plane stress relations: 




T(r) r dr 



T(r)r dr - T (p) 


(B.2) 


0 


rr 




T(r)r dr 


(B.3) 


with the steady-state temperature distribution dlflned In equation (B.l). 


Assuming that the disk Is at a uniform temperature when It Is free of 
stress, and that T(r) from equation (B.l) represents the change In temper- 
ature, aT, substituting equation (B.l) Into equations (B.2) and (B.3) gives: 


and 



(B.4) 


(B.5) 


These stress relations can be Integrated across each disk element to estimate 
an average thermal stress acting on the element. These thermal stresses can 
be superimposed with the disk centrifugal stresses to determine the total ele- 
ment stress condition. Knowing these stresses and stressed volume, the ele- 
mental volume life can be determined. 
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Figure 1.- Analysis methodology flowchart. 
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Figure 2.- Typical finite element model, nominal disk configuration: 
DISK DIAMETER 0.61 M (24 IN.); DISK SPEED, 9000 RPM; NUMBER OF BOLT 
HOLES, 12; BOLT HOLE DIAMTER, 12.7 MM (0.5 IN.); MATERIAL STRESS- 
LIFE EXPONENT, C, 9. 



0 12 3 4 5 

MATERIAL WE I BULL SLOPE 

Figure 3.- Effect of Weibull slope on disk L-jq 

LIFE (DATA NORMALIZED TO: SLOPE = 1. DISK DIAM = 
0.61 M (24.0 IN.), SPEED = 9000 RPM, DISK THICK- 
NESS = 0.0254 M (1 IN.), STRESS LIFE EXPONENT = 
9.0). 
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DISK THICKNESS, IN. 

Figure 5.- Effect of disk thickness on disk L 10 

LIFE (DATA NORMALIZED TO: SLOPE = 1, DISK DIAM 
0.61 M (24.0 IN.), SPEED = 9000 RPM, DISK THICK 
NESS = 0.0254 M (1 IN.), STRESS LIFE EXPONENT = 
9.0). 
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e 
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Figure 6.- effect of speed on disk L 10 life 
(data normalized to: slope = 1. disk dia- 
meter = 0.61 M (24.0 IN.). SPEED = 9000 RPM. 
DISK THICKNESS = 0.0254 M (1 IN.), STRESS 
LIFE EXPONENT = 9.0). 
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MATERIAL STRESS-LIFE EXPONENT 

Figure 7.- Effect of stress-life exponent on 

DISK L 10 LIFE (DATA NORMALIZED TO: SLOPE = 1. 
DISK DIAM = 0.61 M (24.0 IN.), SPEED = 9000 RPM 
DISK THICKNESS = 0.0254 M (1 IN.), STRESS LIFE 
EXPONENT = 9.0). 
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MATERIAL FATIGUE LIMIT, psi 

Figure 8.- Effect of fatigue limit on disk L^q 

LIFE (DATA NORMALIZED TO: SLOPE = 1, DISK DIA- 
METER = 0.61 M (24.0 IN.), SPEED = 9000 RPM, 

DISK THICKNESS = 0.0254 M (1 IN.), STRESS LIFE 
EXPONENT = 9.0). 



Figure 9.- Effect of temperature difference on 

DISK L 10 LIFE (DATA NORMALIZED TO: SLOPE = 1, 
DISK DIAMETER = 0.61 M (24.0 IN.), SPEED = 
9000 RPM, DISK THICKNESS = 0.0254 M (1 IN.), 
STRESS LIFE EXPONENT = 9.0). 







NORMALIZED RELATIVE DISK LIFE. 
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